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Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by various autoantibod- 
ies that recognize autoantigens displayed on the surface of cells undergoing apo ptosis. The genetic con- 
tribution to SLE susceptibility has been widely recognized. We previously reported evidence for linkage 
to SLE of the human chromosome Iq41-q42 region and have now narrowed it from 15 to 5 cM in an 
extended sample using multipoint linkage analysis. Candidate genes within this region include (a) 
PARP, poly(ADP-ribose) polymerase, encoding a zinc-finger DNA-binding protein that is involved in 
DNA repair and apoptosis; (b) TGFB2, encoding a transforming growth factor that regulates cellular 
interactions and responses; and (c) HLX1 } encoding a homeobox protein that may regulate T-cell devel- 
opment. Using a multiallelic, transmission-disequilibrium test (TDT), we found overall skewing of 
transmission of PARP alleles to affected offspring in 124 families (P - 0.00008), preferential transmis- 
sion of a PARP allele to affected offspring (P = 0.0003), and lack of transmission to unaffected offspring 
(P - 0.004). Similar TDT analyses of TGFB2 and HLX1 polymorphisms yielded no evidence for associa- 
tion with SLE. These results suggest that PARP maybe (or is close to) the susceptibility gene within the 
chromosome Iq41-q42 region linked to SLE. 

J. Clin. Invest. 103:1135-1140(1999). 



Introduction 

As is true for many multifactorial complex human dis- 
eases, cumulative studies suggest that interactions of 
multiple genes and environmental factors result in sus- 
ceptibility to systemic lupus erythematosus (SLE) (1,2). 
The prevalence of SLE in the US general population is 
approximately 15-50 per 100,000 (3). The relatively high 
incidence (10-16%) of more than one case in a family has 
suggested a genetic basis for SLE (3). The concordance 
rate of SLE in monozygotic twins (24-57%) is approxi- 
mately 10 times higher than the rate in dizygotic twins 
(2-5%) (1, 3), again implicating genetic factors in this dis- 
order. Based on these epidemiological studies, the rela- 
tive risk for siblings of SLE patients, Xs, defined as the 
risk for siblings of patients divided by the risk in the gen- 
eral population (4), is in the range of 20-40. 

Case-control studies have shown that MHC class II 
genes (HLA-DR and/or DQ but not DP) and certain class 
III genes (C2, C4 } TNFa, and HSP70-2 alleles) confer sus- 
ceptibility to SLE (1, 2). The overall relative risk associ- 
ated with each MHC class II haplotype in various ethnic 
groups is estimated to be two- to threefold. Because of 



extensive linkage disequilibrium within the MHC region, 
the precise loci (and alleles) responsible for disease sus- 
ceptibility have yet to be determined (1, 2). Other non- 
MHC genes have also been associated with SLE. Among 
them, evidence for homozygous deficiency of the first 
complement of the complement system, Clq } predispos- 
ing to SLE is particularly compelling; 90% of such indi- 
viduals have SLE, and Clq knockout mice display an 
SLE-like phenotype (5, 6). In addition, polymorphisms 
in many genes encoding molecules with relevant 
immunological functions (including T-cell receptor a 
and P chains, immunoglobulin allotypes, FcyRIIa, FcyRI- 
Ila, IL-6, IL-10, BCL-2, and mannose-binding protein), 
as well as deletions of specific variable gene segments of 
immunoglobulin genes, have been associated with SLE 
demonstrated mostly by the case-control studies (1, 2, 
7-9). The contribution of these genes to the develop- 
ment of SLE is under investigation. 

Linkage analysis for susceptibility loci of SLE has been 
more extensively studied in murine than in human SLE 
(reviewed in ref. 2). Using the identified murine suscep- 
tibility loci (the overlapping Slel/Nba2/Lbw7) as a guide, 
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we examined seven markers located on a likely syntenic 
human chromosomal Iq31-q42 region in 52 affected 
sibpairs from Caucasian, Asian, and African-American 
families (2, 10-14). Five markers located at Iq41-q42 
showed evidence for linkage (P = 0.0005-0.08) by the 
allele-sharing method using the model-free SIBPAL pro- 
gram of the S.A.G.E. package (ref. 14; Department of 
Biometry and Genetics, Louisiana State University Med- 
ical Center, New Orleans, USA). Subsequently, two inde- 
pendent reports supported linkage of this region; a lod 
score of 1.51 at lq42 using a model-free linkage analysis 
of 105 affected sibpairs, and a lod score of 3.50 at lq41 
using a model-based linkage analysis in 31 African- 
American pedigrees (15, 16). 

In this report, we describe results of fine mapping 
including a multipoint linkage analysis that narrows the 
linked region to 5 cM, and testing three positional can- 
didate genes (PARPj TGFB2, z.ndHLXl) for association 
with SLE. PARP is activated by DNA strand breaks intro- 
duced in the genome through environmental injuries 
(17). It catalyzes the attachment of ADP-ribose units 
from NAD to an array of nuclear acceptor proteins 
involved in cellular proliferation, differentiation and 
DNA repair (17). The potential involvement of PARP in 
SLE has been suggested by its lower-than-normal levels 
of activity and mRNA in patients with SLE and by its 
intermediate range of activity in unaffected family mem- 
bers of patients with SLE (18-20), making it an excellent 
positional candidate. Because of the well-established 
regulatory role ofTGFpl (located on 19ql3) in autoim- 
mune diseases and the fact that same surface receptors 
bind to both TGFpl and TGFP2 (21), we chose TGFB2 as 
the second positional candidate gene. HLX1 resembles a 
diverged Drosophila homeobox gene and is expressed in 
hematopoietic progenitors and activated lymphocytes. 
Functionally, this human homeodomain-containing 
protein has been shown to regulate the development of 
CD4 + T cells in the transgenic mouse model (22). 
Because modulation of CD4 + T cells by administration 
of a specific monoclonal antibody prevents and amelio- 
rates murine SLE (23, 24), we chose HLX1 as the third 
positional candidate gene. Our results of the family- 
based transmission-disequilibrium test (TDT) (25, 26) 
indicate that only one of these three candidate genes, 
PARP, is associated with SLE in this sample. 

Methods 

Subjects. This study was approved by the Human Subject Pro- 
tection Committee of the University of California- Los Angeles 
(UCLA). Multiplex families were recruited by ascertaining 
nuclear families for two or more patients with SLE, which were 
extended to include parents and other siblings. These families 
were recruited at UCLA and at other collaborating sites, includ- 
ing the University of Texas-Houston (F.C. Arnett), the Zen- 
tralkrankenhaus Reinkenheide (K. Hartung), and The Ottawa 
Hospital (R. Goldstein). Simplex families were mainly recruit- 
ed through UCLA rheumatologists to include patients with 
SLE, their parents, and other unaffected sibs if available. Med- 
ical records for patients with SLE were either reviewed by UCLA 
rheumatologists or were provided by the patients' rheumatol- 
ogists in the form of completed checklists (27), and the check- 
list information was then entered into the UCLA database. This 
database was used for classification of patients as SLE (requir- 



ing at least four of the 1 1 American College of Rheumatology 
Classification criteria) (28, 29). Siblings of patients with SLE 
and who have no positive responses in the Connective Tissue 
Screening Questionnaire (30) are classified as unaffected sibs. 

Genotyping. DNA was isolated from peripheral blood cells. 
Microsatellite markers at or near the specific candidate chro- 
mosomal region were selected based on a composite map 
(accessible at http://cedar.genetics.soton.ac.uk/pub/chroml). 
The primers for these markers were purchased from Research 
Genetics Inc. (Huntsville, Alabama, USA). Microsatellite geno- 
typing was deterrnined by scoring the size of PCR products. 
PCR was performed in a 96-well plate using a thermocycler that 
was programmed for 95 °C for one minute, 58 °C for 30 s, and 
72 °C for one minute, with 25 cycles. The labeled PCR products 
were denatured and separated on a 5% sequencing gel. 

The location of candidate genes were based on the GeneMap 
'98 (http://www.ncbi.nlm.nih.gov/genemap) and the Genome 
Database (http://gdbwww.gdb.org). The primer sequences used 
for the determination of genotypes of PARP, HLX1, and TGFB2 
were as described in the Genome Database. The PCR condition 
for PARP and TGFB2 polymorphisms was 30 cycles of 93 °C for 
one minute, 56° C for 30 s, and 72 °C for one minute, plus a 
final extension of 72 °C for five minutes. The HLX1 PCR con- 
dition was 10 cycles of94°Cfor 15 s,55°C for 15 s, and 72°C 
for 30 s, followed by 20 cycles of 89° C for 15 s, 55° C for 15 s, 
and 72 °C for 30 s, with a final extension at 72 °C for 10 min. 
For these two candidate genes, the PCR mixture contained 40 
ng genomic DNA, 20 mM Tris-HCl (pH 8.0), 2 mM MgCl 2 , 10 
mM KC1, 6 mM (NH 4 ) 2 S0 4 , 0.1% Triton X-100, 10 Hg/ml BSA, 
0,5 U native Pfu DNA polymerase (Stratagene, San Diego, Cal- 
ifornia, USA), 200 |i.M dNTP, and 0.2 uM primers in a 5-jll reac- 
tion. The PCR mixture for TGFB2 was similar to that described 
for PARP and HLX1, except for the buffer, which contained 10 
mM Tris-HCl (pH 9.0), 1.5 mM MgCl 2 , 50 mM KC1, and 0.1% 
Triton X-100. Aliquots of fluorescent PCR products of PARP 
and HLX1 were electrophoresed using a 377 Prism ABI 
sequencer and analyzed by GeneScan and Genotyper programs 
(Applied Biosystems, Foster City, California, USA). The TGFB2 
genotyping was determined using radioactively labeled primers 
and autoradiography. 

Data analysis. Multipoint linkage analysis of 14 microsatellite 
markers in the chromosome lq31-42 region and the SLE phe- 
notype (affected or unaffected) was conducted on nuclear fam- 
ilies using the MAPMAKER/SIBS program (31) on a Sun Ultra 
60 workstation (Sun Microsystems Inc., Palo Alto, California, 
USA), with a UNIX operating system. The program simultane- 
ously incorporates allele-sharing information at multiple mark- 
ers in a chromosomal region to estimate the expected propor- 
tions of affected sibpairs sharing zero, one, and two alleles 
identical by descent for points along the chromosomal region. 
Spouses of probands of multiplex families were used to esti- 
mate allele frequency. A lod score is generated by taking the log 
of the likelihood of the data at the maximum likelihood esti- 
mate of allele-sharing proportion, compared with the likeli- 
hood under the expectation of no linkage to the region. 

The TDT can detect an allelic association of a marker with 
disease in the presence of its linkage to the disease (25). The 
method involves assessing the transmission of alleles from het- 
erozygous parents to their affected offspring and testing 
whether there is evidence that transmission of alleles at a mark- 
er is different from 0.5 for each allele. To assess the deviation 
from 0.5 for all alleles simultaneously, the extended TDT 
(ETDT), which uses a logistic regression approach, has been 
conducted (software provided by the program authors) (26). If 
a significant deviation from 0.5 was detected, the same software 
was employed to examine each allele separately for a deviation 
from a 0.5 probability of transmission. 

The relative risk for the 85-bp allele identified by the ETDT at 
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Figure 1 

Multipoint linkage analysis of the chromosome 
1q31-q42 region using the MAPMAKER/SIBS pro- 
gram. Positions of genetic markers relative to 
D1S510 are expressed in cM as the relative genet- 
ic distance shown on the x axis. All markers except 
HLX1 and TCFB2 were used in this analysis. The 
heuristic guideline of 1 lod below the peak value 
was used to identify a confidence interval of 5 cM 
(flanked by D1S2860 and D1S213) for the loca- 
tion of the putative SLE susceptibility gene. 
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the PARP locus has been assessed using a logistic regression 
method (32). Poisson regression analysis was applied to the num- 
bers of parent-child triads observed among the oldest affected off- 
spring with at least one parent heterozygous for the 85-bp allele. 
The risk of developing SLE with one or two 85-bp alleles was esti- 
mated by using the GENMOD procedure of the Statistical Analy- 
sis System (SAS) package computer program (33). 

Results 

We previously reported evidence for linkage of a 15-cM 
region on chromosome 1 to SLE using the affected sib- 
pair approach (14). To support and extend this finding, 
we have increased our sample size by 50%, from the ini- 
tial 52 to 78 SLE-affected sibpairs from Caucasian, 
Asian, and African-American families. All participating 
members (77 parents; 134 affected and 89 unaffected 
siblings) were genotyped using 14 microsatellite mark- 
ers within a 30-cM region (Figure 1) containing the pre- 
viously identified 15-cM region. Previously, this linked 
region was established by model-free two-point linkage 
analyses using seven markers (D1S510, D1S249, 
D1S245, D1S229, D1S213, D1S225, and D1S103). As 
shown in Figure 1, multipoint linkage analysis on the 
current sample using the model-free MAPMAKER/SIBS 
program (31) identified a peak with a lod score of 3.3 
and narrowed the region of interest to approximately 5 
cM between markers D1S2860 and D1S213. 

Because positive association results of candidate genes 
located within identified linked regions can narrow the 
search for susceptibility loci (e.g., apolipoprotein E type 
4 allele in late-onset familial Alzheimer's disease and 
CTLA4 in type 1 diabetes), we tested three candidate 
genes with known polymorphisms that are likely to 



localize within this narrowed region — specifically PARP 
(also known as ADPRT), TGFB2, and HLX1 - using a 
multiallelic TDT approach (25, 26). The PARP polymor- 
phism we tested is a CA repeat located 906 bp upstream 
of the transcription start site (36). There were nine PARP 
alleles detected in our cohort of 124 multiplex and sim- 
plex families consisting of 75 Caucasian, 25 Hispanic, 17 
Asian, and seven African-American families, containing 
223 parents, 289 patients, and 226 unaffected sibs. The 
oldest affected and unaffected offspring in each family 
were used for these analyses, as other sibs from the same 
families are not independent samples. Because racial ori- 
gin is likely to contribute to genetic heterogeneity of 
SLE, we compared these candidate genes using one eth- 
nic group, our largest one: the Caucasian sample. Table 
1 summarizes data from all PARP alleles analyzed by the 
multiallelic TDT (25, 26) using our multiethnic cohort, 
and Table 2 summarizes data from the 75 Caucasian 
families. Analysis of the combined four frequent alleles 
(85, 93, 95, and 97 bp) showed the skewed transmission 
to affected offspring in the combined sample containing 
all ethnic groups (P = 0.00008; Table 1) as well as in the 
Caucasian sample (P - 0.002; Table 2). Further inclusion 
of the five infrequent PARP alleles (87, 89, 91, 99, and 
101 bp) in the multiallelic TDT yielded an overall P value 
of 0.001, supporting significantly disproportional trans- 
mission to affected offspring when all nine alleles were 
analyzed as a group. Other ethnic groups were each too 
small in number to be analyzed separately. 

Application of the TDT to these PARP alleles indicat- 
ed that the 85-bp allele of PARP was preferentially trans- 
mitted to affected offspring (transmission/nontrans- 
mission = 66:33; P = 0.0003) and preferentially not 
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Table 1 

Skewed transmission of PARP alleles in the multiethnic cohort 

All ethnic groups combined 



Alleles 




Affected offspring 






Unaffected offspring 


(bp) 


t 


nt 


%t 


P value 


t 


nt 


%t 


P value 


85 


66 


33 


67 


0.0003 


24 


52 


31 


0.004 


93 


28 


36 


44 


NS 


24 


21 


53 


NS 


95 


13 


15 


46 


NS 


13 


5 


72 


NS 


97 


19 


42 


31 


0.0009 


31 


20 


61 


NS 


87 


4 


5 






5 


3 






89 


0 


3 






2 


0 






91 


1 


1 






0 


1 






99 


4 


0 






2 


0 






101 


1 


1 






1 


0 






Multialleles 






0.00008 








0.02 



Data presented are derived from all PARP alleles transmitted (t) and not transmitted (nt) from heterozygous 
parents to offspring in 1 24 families. Data are shown as the percentage transmission of each allele (%t). The 
P value for each allele is the level of significance based on the deviation from the expected random (50%) allele 
transmission and was evaluated using a chi-squared statistic (x 2 - [t - nt) 2 /(t + nt)] with 1 degree of freedom). 
The multiallelic P value is the level of significance based on testing the four frequent alleles (85, 93, 95, and 97 
bp) simultaneously for a deviation from the expected 50% transmission. NS, not significant (P > 0.05). Indi- 
vidual P values for the infrequent alleles (87, 89, 91 , 99, and 1 01 bp) were not calculated. 



Table 2 

Skewed transmission of PARP alleles in Caucasians 



showing no obvious distortion from the 
expected random distribution (% 2 = 
0.36; P = 0.6). TDT analyses of HLX1 
were also performed using an intronic 
dinucleotide repeat within this gene in 
an analogous fashion, but they yielded 
no evidence for association with SLE. 
Transmission data of the seven HLX1 
alleles to affected offspring from 129 
heterozygous Caucasian parents are 
summarized in Table 3; no significant 
difference from the expected random 
distribution was detected (% 2 for the 
genotype-wise TDT = 1 1.35; P = 0.88). In 
summary, TDT results of these three 
candidate genes yielded only evidence 
for association of the PARP polymor- 
phism with SLE in our tested cohort. 
Thus, our data suggest that PARP may 
be the SLE susceptibility gene within 
the linked region previously identified, 
or at least is close to it. 



Caucasian 



Alleles 




Affected offspring 






Unaffected offspring 




(bp) 


t 


nt 


%t 


P value 


t 


nt 


%t 


P value 


85 


41 


17 


71 


0.001 


16 


32 


32 


0.03 


93 


14 


18 


44 


NS 


12 


9 


57 


NS 


95 


6 


11 


35 


NS 


10 


2 


83 


0.03 


97 


10 


23 


30 


0.02 


16 


12 


57 


NS 


87 


1 


2 






2 


1 






89 


0 


0 






0 


0 






91 


0 


1 






0 


1 






99 


0 


0 






0 


0 






101 


1 


1 






1 


0 






multialleles 








0.002 








NS 



Data presented are derived from 75 Caucasian families. See Table 1 for definitions oft, nt, %t, NS, and P value. 



transmitted to unaffected offspring (transmission/non- 
transmissiow = 24:52; P = 0.004) in the 124 multiethnic 
families (Table 1). These data also indicated that the 97- 
bp allele of PARP might be protective, as it was preferen- 
tially not transmitted to affected offspring (transmis- 
sion/nontransmissiow = 19:42; P = 0.0009) (Table 1). 
Consistently, in Caucasian families (Table 2), the 85-bp 
allele of PARP was significantly transmitted to offspring 
with SLE (P = 0.001), whereas the 97-bp allele was trans- 
mitted more frequently to unaffected offspring (P = 
0.02). Compared with having no 85-bp alleles, the risk of 
developing SLE with one is 2.3 and is 4.0 for being 
homozygous (32). 

Similar TDT analyses of a dinucleotide repeat within 
the TGFB2 was performed using 81 Caucasian nuclear 
families. However, transmission data of seven TGFB2 
alleles from 149 parents to affected offspring demon- 
strated homozygosity in 83% of the analyzed parents. 
Among the 25 heterozygous parents studied, the most 
frequent TGFB2 allele was transmitted 11 times, but 
not transmitted 14 times, to affected offspring, thus 



Discussion 

In this study, we have demonstrated that 
the 85-bp PARP allele is preferentially 
transmitted from heterozygous parents 
to offspring affected with SLE in a sam- 
ple of 124 families comprising Cau- 
casians, Hispanics, Asians, and African- 
Americans. Our rationale for performing 
the TDT analysis is to identify a genetic 
marker at or near the putative suscepti- 
bility gene within the large candidate 
region suggested by linkage analyses. 
Evidence for linkage of chromosome 
Iq41-q42 in human SLE is very strong, 
including (a) our previous report of a 15- 
cM region established with 52 affected 
sibpairs from multiple ethnic groups 
(14); (b) the narrowed 5-cM region with 
a maximum lod score of 3.3, as shown in Figure 1; (c) the 
lod score of 1.51 of a marker at lq42 in an independent 
sample of 105 affected sibpairs (15); and (d) the lod score 
>2 at lq41 in 31 African- 
American pedigrees (16). 
In the presence of estab- 
lished linkage to this 
region, we have tested 
three positional candidate 
genes for association by 
using the oldest affected 
offspring of each nuclear 
family in the TDT analy- 
sis. This family-based asso- 
ciation test is considered 
to be a better study design 
than the population-based 
association test because 
the latter is more prone to 
spurious associations if 
cases and controls are 



Table 3 

HLX1 alleles are not associated with 
susceptibility to SLE 



Alleles 




Affected offspring 




(bp) 


t 


nt 


%t 


168 


9 


8 


53% 


166 


10 


9 


53% 


164 


8 


11 


42% 


162 


17 


18 


49% 


160 


30 


30 


50% 


158 


39 


32 


55% 


152 


16 


21 


43% 



The analyzed Caucasian families contained 
129 heterozygous parents. The overall 
transmission of 7 tested alleles was not dif- 
ferent from the expected random distribu- 
tion with a P value of 0.88. See Table 1 for 
definitions of t, nt, %t, NS and P value. 
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drawn from unmatched populations (37). That none of 
these three candidate genes are located at the peak of the 
observed linkage (Figure 1) does not negate the possibil- 
ity of their being the susceptibility gene. This may main- 
ly reflect the imprecision of fine mapping using the link- 
age method due to the lack of one-to-one correlation 
between genotypes and phenotypes in complex diseases 
(37). These tested positional candidates have known 
polymorphisms of dinucleotide repeats with maximum 
heterozygosity indexes of 0.63, 0.52, and 0.82 for PARP y 
TGFB2, and HLX1, respectively (38). The low heterozy- 
gosity index of this polymorphism within TGFB2 
explains the low percentage of heterozygous parents we 
observed. Our studies of these three candidates showed 
that only PARP alleles were associated with SLE. 

Our positive TDT results suggest that we are very close 
to the susceptibility gene based on linkage disequilibri- 
um. The extent of linkage disequilibrium in the human 
genome appears quite variable; from 1-2 cM as observed 
in the MHC region to within 4.5 kb from studies of 
markers in the insulin gene (1, 39). In general, linkage 
disequilibrium may exist within a few hundred kilobase 
in most human genomic regions of nonisolated popula- 
tions (40). Because map positions of PARP, TGFB2, HLX1 
and their flanking genes are all currently estimated to be 
within 6 megabase (38, 41), detailed physical mapping of 
the region surrounding PARP will be necessary to pro- 
vide additional positional candidate genes. 

An equally likely interpretation of our TDT results is 
that PARP is the SLE susceptibility gene within this 
linked region. The polymorphism of PARP we tested is 
located in the promoter region near the binding site of 
the transcription factor Ying Yang 1 (42) and, thus, may 
affect transcription. The functional role of PARP has 
been extensively studied by using a specific inhibitor (3- 
aminobenzamide or antisense RNA) and by studies of 
knockout mice (43). Cumulative data have shown that 
the absence of PARP activity results in elevated sponta- 
neous genetic rearrangements and hypersensitive 
responses to DNA damage, suggesting a role for PARP in 
maintaining genomic stability (43). In addition, PARP is 
known to be involved in induced apoptosis. Specific 
cleavage of PARP occurs at an early stage of apoptosis, 
which may commit cells for self-destruction by elimi- 
nating protective mechanisms, including DNA repair 
(17). Both defective DNA repair and abnormal apopto- 
sis have been implicated as causative factors for SLE 
(44-46). Although no gross defects in apoptosis are 
found in PARP knockout mice, splenocytes of these mice 
display a more rapid apoptotic response to an alkylating 
agent (47). Cell lines with disrupted PARP expression 
show insensitivity to apoptotic signals (48). Both exper- 
imental systems suggest a regulatory role of PARP in 
induced apoptosis, with impaired apoptosis less 
detectable in whole animals than in cell lines, probably 
because of other compensatory routes. Unlike Clq 
knockout mice, lupus-like phenotypes have not been 
reported in PARP knockout mice. It is worth noting that 
only 25% of Clq-deficient mice in the appropriate genet- 
ic background have histological evidence of glomeru- 
lonephritis at eight months of age (6). This observation 
has not been formally tested in PARP-deficient mice at 



the appropriate age. Our working hypothesis is that the 
85-bp allele of PARP (or polymorphisms in linkage dise- 
quilibrium within PARP) confers defective DNA repair 
and abnormal apoptosis and thus predisposes to SLE. 
Previous studies of PARP activity and mRNA in patients 
with SLE and in family members further suggest a direct 
role of PARP in susceptibility for SLE (18-20). Our link- 
age and association suggest that PARP is an SLE suscep- 
tibility gene or, at the very least, is in physical proximity 
to an SLE susceptibility gene. 
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